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The electromagnetic two-body problem has neutral differential delay equations of motion that, 
for generic boundary data, can have solutions with discontinuous derivatives. If one wants to use 
these neutral differential delay equations with arbitrary boundary data, solutions with discontinu- 
ous derivatives must be expected and allowed. Surprisingly, Wheeler-Feynman electrodynamics has 
a boundary value variational method for which minimizer trajectories with discontinuous deriva- 
tives are also expected, as we show here. The variational method defines continuous trajectories 
with piecewise defined velocities and accelerations, and electromagnetic fields defined by the Euler- 
Lagrange equations on trajectory points. Here we use the piecewise defined minimizers with the 
Lienard-Wierchert formulas to define generalized electromagnetic fields almost everywhere (but on 
sets of points of zero measure where the advanced/retarded velocities and/or accelerations are dis- 
continuous). Along with this generalization we formulate the generalized absorber hypothesis that 
the far fields vanish asymptotically almost everywhere and show that localized orbits with far fields 
vanishing almost everywhere must have discontinuous velocities on sewing chains of breaking points. 
We give the general solution for localized orbits with vanishing far fields by solving a (linear) neutral 
differential delay equation for these far fields. We discuss the physics of orbits with discontinuous 
derivatives stressing the differences to the variational methods of classical mechanics and the exis- 
tence of a spinorial four-current associated with the generalized variational electrodynamics. 
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INTRODUCTION 

Non-radiating motion of extended charge distributions 
in classical electrodynamics has been known to exist for 
some time (c.f. [ll-Q and references therein, and (a-Q)- 
On the other hand, for systems with a few point charges, 
Larmor's radiation of energy at a rate proportional to 
the squared modulus of the acceleration plagues classical 
electrodynamics. To construct orbits that do not radi- 
ate, and hence are without acceleration, a simple option 
are constant velocity motions, which imply unbounded 
motion. 

Along bounded two body motions supported by mu- 
tual action at a distance, we expect acceleration to be 
needed to change velocities, unless velocities are allowed 
to change discontinuously. For example, periodic polygo- 
nal orbits with piecewise constant velocity segments have 
vanishing radiation fields. 

Here we extend Wheeler-Feynman electrodynamics [8[ 
to include motion with discontinuous velocities. This is 
a natural extension provided by the variational bound- 
ary value problem [9\. The resulting extended electro- 
dynamics has several appealing physical features: (i) 
There exists a scalar function (the finite action [9|), and 
the condition for a minimizer demands that the partial 
derivatives of the action, with respect to each particle's 
four-velocity, be continuous along minimal orbits. These 
continuous four-component linear currents are analogous 
to the Dirac-equation of quantum mechanics, thus en- 



dowing the extended Wheeler-Feynman electrodynamics 
with spin. This is a feature not present in any other 
classical electrodynamics of point-charges; (ii) Besides 
naturally including non-radiating orbits, the extended 
electrodynamics can be shown to lead simply to a de 
Broglie length for double-slit scattering upon detailed 
modeling [10|; (iii) The absorber hypothesis, first ide- 
alized to hold as an average over an infinite universe [8|, 
has no known solutions 11 ill for many-body motion in 
Wheeler-Feynman theory |llMl5| with which it is consis- 
tent. Here we show that the variational electrodynamics 
allows a concrete realization of the absorber hypothesis 
for a two-particle universe, i.e., there exists a non-empty 
class of two-body motions with vanishing far-fields, so 
that we do not need either large universes or random- 
ization [la . llTJ I ; and (iv) two-body orbits with vanishing 



far-fields were used in Ref. [18| to predict spectroscopic 



lines for hydrogen with a few percent precision. 

Since the speed of light is constant in incrtial frames, 
the equations of motion for point-charges are state de- 
pendent differential delay equations. More specifically, 
Wheeler-Feynman electrodynamics [8|, [lOl, llSJ has mixed- 
type state-dependent neutral differential delay equations 
of motion for the two-body problem. 

The theory of delay equations is still incomplete [20, 



2l[ but it is known that purely-retarded differential delay 



equations with generic C^ initial histories have continu- 
ous solutions with a discontinuous derivative at the ini- 
tial time. The derivative becomes continuous at the next 
breaking point [20| and progresses from C'^ to C'^+^at 



successive breaking points. On the other hand, a purely 
retarded neutral differential delay equation with a generic 
C^ initial history [20j can have continuous solutions with 
discontinuous derivatives at all breaking points. 

If one wants to use the electromagnetic neutral differ- 
ential delay equations with arbitrary boundary data, so- 
lutions with discontinuous derivatives must be expected 
and accommodated. Surprisingly, this same neutrality is 
compatible with the recently developed boundary-value- 
variational method for Wheeler- Fey nman electrodynam- 
ics [9| . For orbits where the acceleration is not defined at 
a few points, the variational method offers a well-posed 
alternative to define trajectories beyond those satisfying 
a Newtonian-like neutral differential delay equation ev- 
erywhere. The variational method involves an integral 
that requires only piecewise-defined velocities, general- 
izing naturally to continuous orbits with discontinuous 
derivatives at breaking points. 

Our generalized electrodynamics contains the C^ or- 
bits of the Wheeler- Fey nman theory. As shown in Ref. 
y, if boundary data are such that the extremum orbit 
is piecewise C^ with continuous velocities, the Wheeler- 
Feynman equations hold everywhere with the exception 
of a countable set of points where accelerations are dis- 
continuous (which is a set of measure zero for the ac- 
tion integral). We henceforth define a breaking point 
as a point where velocity or acceleration are discontinu- 
ous. Here we show that continuous orbits with discon- 
tinuous velocities are possible minimizers if these sat- 
isfy further continuity-conditions. These continuity con- 
ditions are non-local, unlike the conditions for an ex- 
tremum of the variational methods of classical mechan- 
ics, which do not allow discontinuous velocities. Fi- 
nally, if the extremum is not piecewise C^, the varia- 
tional method defines minimizers that are not described 
by piecewise-defined- Wheeler-Feynman neutral differen- 
tial delay equations (which are not studied here). 

To discuss the relationship to Maxwell's electrodynam- 
ics it is important to keep in mind that: (i) Wheeler- 
Feynman electrodynamics is a theory of trajectories, 
where fields are only derived quantities; and (ii) the 
boundary-value-variational-method defines only a finite 
segment of a trajectory, rather than a global trajectory 
[9|. The variational equations along piecewise C^ orbits 
include the electromagnetic fields in the Euler-Lagrange 
equations |9|], which are used here to give a derived op- 
erational meaning to the electromagnetic fields [22]. The 
electromagnetic fields appear as coupling terms of the 
variational equations and are defined on trajectory seg- 
ments by the usual electromagnetic formulae [9|. 

In our generalization we use the Lienard-Wierchert 
electromagnetic formulae to define fields by extension at 
all space-time-points for which future and past lightcones 
fall in the finite segment of the minimizer trajectory. 
For continuous trajectories with discontinuous velocities, 
and/or accelerations on sets of measure zero, we con- 



struct the electromagnetic fields only for points having 
a future and past lightcone, leaving the fields undefined 
where the past or future lightcones have a discontinu- 
ous velocity/acceleration (usually another set of measure 
zero). We further introduce the concept of short-range 
orbits as localized orbits with far-fields vanishing almost 
everywhere. This bears a close relation to the electro- 
magnetic notion of radiation [22]. 

In their original articles, Wheeler and Feynman [8J at- 
tempted to derive an electrodynamics with retarded-only 
fields from the hypothesis that the universal far-fields 
vanish at all times (the absorber hypothesis) ^. Here 
we generalize the absorber hypothesis [g to include fields 
that can be undefined on sets of measure zero, thus ar- 
riving at the generalized absorber hypothesis (GAH) that 
the far-fields vanish almost everywhere. We show that 
short-range-two-body-orbits must involve discontinuous 
derivatives on a countable set of points. 

One advantage of our generalization is to include 
spatially-bounded-globally-defined-continuous-orbits 
with far-fields vanishing almost everywhere, which we 
call short-range orbits. This generalization presents 
itself naturally as the next option after one shows that 
there are no C^ localized orbits with far-fields vanishing 
everywhere. The short-range piecewise C^ continuous 
orbits are naturally described by the variational method 
in the same way as the globally C^ continuous orbits. 
However, the former are minimizers inside a larger 
family of boundary-data, which is the second advantage 
of our generalization. This extended class of orbits 
includes orbits that are limits of Cauchy sequences of 
orbits with far-fields disturbing the universe less and 
less, i.e., the vanishing-far-field-limit of the GAH. 

In this paper we use the word "minimizer" meaning a 
generalized critical point of the variational method, that 
could be either a minimum or a saddle point. The paper 
is divided as follows: In Section II we discuss the vari- 
ational method for piecewise-defined continuous orbits 
with discontinuous derivatives. We show that the vari- 
ational method prescribes a continuous momentum cur- 
rent at each breaking point in addition to Euler-Lagrange 
equations from each side of the breaking point. We dis- 
cuss how the non-local-momentum-currents can be con- 
served even in the presence of velocity discontinuities 
along "sewing chains" . In Section HI we prove that 
globally-defined short-range bounded orbits must have 
discontinuous velocities on a sewing chain of breaking 
points by giving the general solution to a neutral differ- 
ential delay equation for the far-fields. In Section IV we 
discuss the physics of generalized minimizers along with 
some open questions and differences from bounded orbits 
to unbounded scattering orbits. 



BOUNDARY VALUE VARIATIONAL METHOD 



The variational method [9| is well defined for continu- 
ous trajectories Xi{ti) and X2{t2) £ R^ that are piecewise 
C^ . The boundary conditions for the variational method 
(9| are illustrated in Figure 1, i.e., the initial point Oa for 
trajectory 1 plus the segment of trajectory 2 inside the 
lightcone of Oa, and the endpoint Lb for the trajectory 
2 plus the segment of trajectory of particle 1 inside the 
lightcone of Lb ■ For variations of trajectory 1 the action 
functional y\ reduces to 
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where K2 depends only on trajectory 2 and quantities 
of particle 2 are defined at times i2±(xi(ti)) according 
to the implicit condition for the advanced/retarded light 
cones of ti, i.e., 



i2±(xi(ti)) ^h± |xi(ti) - X2(i2±) 



(2) 



In Eq. dl]) the ri2+ = |xi(ii) — X2(i2±)| are the distances 
from xi(ii) to the respective advanced/retarded posi- 
tion X2(i2±) along trajectory 2, unit vector ni2± points 
from Xi(ii) to the respective advanced/retarded position 
X2(t2±), i-e., ni2± = (xi(ti) - X2(i2±))Ai2± and last 
V2± = dyi2/dt2\t2±- Notice that Eq. ([T]) is an inte- 
gral over the velocities, and is a well-defined operation 
even for trajectories with discontinuous velocities (and 
even more general types of continuous trajectories with 
square-integrable velocities that are not studied here) . 

Here we extend Wheeler- Fey nman electrodynamics to 
trajectories with discontinuous velocities on a count- 
able set of points using the boundary- value- variational- 
method. For piecewise C^ trajectories (and piecewise 
C^ histories) it is possible to define disjoint intervals 
t e [l'^_i,l~], with l~ = Z+, for a = l,...,iV, where 
the continuous trajectory Xi{ti) and delayed arguments 
t2±{ti) are piecewise C^. The upper plus in /+ indicates 
the right-limit of the a*'* breaking point while the upper 
minus in l~ indicates the left-limit of the a*^^ breaking 
point. These are not to be confused with the lower in- 
dices used to denote quantities evolving in future and 
past light cones. 

The variations of trajectory 1 for the action ([1]) are 
defined piecewise C^ with fixed endpoints, i.e.. 
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FIG. 1: Illustrated in red is the initial point Oa of trajectory 
1 plus the segment of trajectory 2 inside the lightcone of Oa, 
i.e., from point 0~ to point O^ and the endpoint Lb of 
trajectory 2 plus the segment of trajectory of particle 1 inside 
the lightcone of Lb, i.e., from point L~ to point L"*". The 
trajectory of particle 1 of the variational method goes from 
Oa to L~ (blue line) while the trajectory of particle 2 goes 
from O"*" to Lb (green line). The first breaking point is point 
O+which generates a forward sewing chain of breaking points 
/i,/2,/3 while endpoint L" is a breaking point generating a 
backward sewing chain of breaking points 61 , 62 , fes . 



where and overdot denotes a time-derivative and the 
boundary conditions are 



bi(;+ = o) = o. 



(4) 



If the continuous and piecewise C^ perturbation bi(ii) 
has a discontinuous derivative in another set of intervals 
t S {ht-iihT^), then the perturbed trajectory Ui(ti) is 
continuous and piecewise C^ in the extended set of inter- 
vals defined by all intersections of the sets {h't_i,h~) and 
(^ct_i,^ct)- This simply increases the number of piecewise 
intervals {l'^_i, l'^) up to a ~ M > N and the boundary 
condition for hi{Tj^-) of Eq. ^ reads 

h^{ll,=TL-) = 0. (5) 

Substituting the perturbed trajectory ([3]) into the ac- 
tion ID) and making a linear expansion about the orbit 
defines the Frechet derivative, i.e.. 
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'9vi 



bi)]dti-fo(| bil), (6) 



where a "•" indicates the scalar product in M'' and | bi | is 
the sup-norm for the Banach space of piecewise C^ varia- 
tions [2i]. In particular if the orbit a;i(ti) : [0,r/^-] — > M^ 
is piecewise C^ then Ui(ii) is continuous and piecewise 
C^ on the same extended set of intervals. 

We can integrate Eq. © by parts in each interval 



yielding 
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(7) 






Since hi(ti) is continuous we can re-arrange the second 
term of the right-hand-side of Eq. ([7]) to give 
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where 






(9) 



Equation (j9]) defines the momentum jump at i = /^T, i.e., 
the first (second) term on the right-hand-side of Eq. ^ is 
the momentum evaluated from the right (left) of i = /~ . 
The conditions for a critical point in the class of con- 
tinuous piecewise C^ orbital variations of the piecewise 
C^ continuous orbit are: (i) satisfy the Euler-Lagrange 
equations piecewise, to make the first term on the right- 
hand-side of Eq. ([5]) vanish; and (ii) have a continuous 
momentum d^/dvi at the breaking points so each term 
of the sum of the right-hand-side of Eq. ([8|) vanishes for 
arbitrary hi{l~) , i.e.. 



"^-(O = l^('i 



(10) 
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2ri2-(l - ni2- • V2-) 2ri2+(l + ni2+ • V2+) 

which displays a surprising difference compared with the 
result obtained from variational methods in classical me- 
chanics. 



As illustrated in Figure 1, a simple piecewise-defined 
orbit has the "sewing chain" of breaking points (/i,/2, . . . 
and &i, 627 • ■ •)> where one velocity can jump if the other 
velocity has jumped at either the past or future break- 
ing point. The two-body-Noether-momentum (formula 
(A23) of Ref. [9|), involves an integral that is insensitive 
to velocity jumps plus two non-local momentum terms 
given by Eq. ^ (see Eqs. (A25) and (A26) of Ref. [9,]), 
that are sensitive to jumps. Therefore, the two-body- 
Noether-momentum is conserved as long as Eq. ([Tl]) is 
continuous across the jumps. The first term on the right- 
hand-side of Eq. ([8|) is the Wheeler- Fey nman equation of 
motion for particle 1 (i,-e., the usual Euler-Lagrange 
equation restricted here to piecewise segments). Mini- 
mization respect to variations of trajectory 2 yields the 
neutral differential delay equation of motion for particle 
2, and an analogous continuity with indices 1 and 2 ex- 
changed. Notice that this surprising difference compared 
with the results from variational principles in classical 
mechanics requires a minimum of two bodies and a non- 
local Lagrangian. 

We note that the electromagnetic variational method 
has a parametrization-invariance-symmetry that allows 
the action ((TJ to be expressed in Minkowski-four-space 
using four-velocities respect to an arbitrary evolution 
parameter |9j]. The derivative of the parametrization- 
invariant-Lagrangian with respect to the first-component 
of the four-velocity is. 



mi 



yr^^ 



(12) 



As is usual in the neighborhood of breaking points one 
defines derivatives from the left-hand-side and from the 
right-hand-side [20|. For the local Lagrangians of clas- 
sical mechanics one usually has dJif/dvi = G'i(xi,vi), 
which combined with Eq. (jlOp along a continuous trajec- 
tory would imply that each velocity is continuous. Con- 
tinuity of velocity along a piecewise C^ continuous orbit 
combined with the Euler-Lagrange equations from each 
side of the breaking point further determine a continu- 
ous acceleration, so that the orbit is actually C^ at the 
breaking point. Therefore, in classical mechanics the re- 
striction to piecewise C^ orbits implies globally C^ or- 
bits. 

However, for Eq. ([T|), or the Lagrangian given in Eq. 
(14) of Ref [9], the continuous momentum term is 

d^ TOiVl 



2ri2-(l - ni2- ■ V2-) 2ri2+(l + ni2+ • V2+) 



'dvf 



Equation (J12l) represents the time-component of the four- 
momentum which must be continuous at the breaking 
points of a minimizing trajectory. This is a generalization 
of the argument leading to Eq. pT|) . There is a four- 
current associated with the minimization respect to each 
particle's trajectory. 

Last, after solving for velocity discontinuities along 
sewing chains of breaking points, the second condition 
for a minimizer are the piecewise-restricted Wheeler- 
Feynman equations of motion. These hold at each side 
of a breaking point and involve the limiting accelerations 
from the two sides of that breaking point. The discon- 
tinuous velocities satisfying Eqs. pTt and P^ . when 
substituted into the Wheeler- Fey nman equations at each 
side of a breaking point, then define a condition to be 
satisfied by the acceleration discontinuities. Since this 
condition involves a singular matrix [30| , acceleration dis- 
continuities are not fully determined by velocity discon- 
tinuities. In general, velocity discontinuities cause accel- 
eration discontinuities, even though there can be special 
orbits with continuous velocities and discontinuous ac- 
celerations along the null direction of the singular ma- 



trix [30|. Only in that case our generalization is equiva- 
lent to piecewise restricted Wheeler-Feynman equations, 
otherwise completion by a finite action [9| includes other 
types of trajectories. In either case, along a piecewise 
defined orbit, continuity of the spinor currents (jlip and 
(|12p ensure the well-posed continuation of the minimizer 
across each breaking point. For example, solutions with 
continuous velocities yield trivially continuous momenta 
(fTTj) and (fT2|) . so that any finite portion of a global solu- 
tion of the piecewise-restricted Wheeler-Feynman equa- 
tions with continuous velocities is a minimizer of the fi- 
nite variational method with suitable boundaries. For 
solutions with discontinuous velocities, continuation to 
a global trajectory is possible using the (discontinuous) 
velocity determined by solving conditions (fTTj) and (fT2|) 
for the most advanced velocity (to be shown elsewhere). 



SHORT-RANGE ORBITS 



Wheeler-Feynman electrodynamics is a theory of direct 
interaction between charges [8|, \2^ . The boundary- value- 
variational-method (previous section and Fig. 1) defines 
minimizers with a vanishing Frechet derivative ^ be- 
tween the time-spans of Fig. 1, rather than globally- 
defined trajectories. As shown in Ref. [9|, the two- 
body-Euler-Lagrange equations can be cast in the form 
of Newtonian equations of motion with each acceleration 
multiplied by the mass on the left-hand-side, while the 
right-hand-side has the form of a Lorentz -force-law. It 
is precisely these Euler-Lagrange-equations that define 
the electromagnetic fields of Wheeler-Feynman theory as 
derived quantities evaluated on trajectories. 

Extending these fields defined by the Lorentz-sector-of- 
Euler-Lagrange-equations to fields on positions outside 
trajectories is tricky, because in a theory of trajectories 
one should: (i) add a third particle to the variational 
problem; and (ii) arrange things such that the third tra- 
jectory passes by the desired point. Obviously, a third 
charge changes the minimization problem and perturbs 
the original two-body-orbit, unless it can be placed so 
far that its couplings to the original two-body-orbit are 
small. A bounded GAH two-body-orbit is special be- 
cause its far-fields vanish almost everywhere and a third 
trajectory can be placed reasonably near without dis- 
turbing the two-body-orbit. Keeping in mind that the 
far-fields are the strongest couplings to a third charge, we 
investigate the existence of such (localized) short-range- 
two-body-orbits (GAH). 

We now adopt a unit system in which the speed of 
light is c = 1 and apply the usual formulae of electrody- 
namics to piecewise-defined-trajectories with the excep- 
tion of points where past/future velocities/accelerations 
are undefined, i.e., the fields are undefined on a set of 
measure zero. We consider continuous piecewise C^ tra- 



jectories Xfc(tfe) enclosed by a sphere of radius R in an 
inertial frame. We specify space-time points (i, Rn) on 
the sphere by a time t and unit vector n normal to the 
surface of the sphere, and introduce an index fc = 1, 2 to 
label the charges. 

The far-electric field of a point charge in the Wheeler- 
Feynman electrodynamics is the sum of the half- 
advanced/half- retarded fields [8|, 



1 



E{t,Rn)^--E 



adv 



-£'■ 



while the far-magnetic field is given by 



B{t,Rn) = -n+ x E 



adv 



1 



-n_ X E"^ 



(13) 



(14) 



The unit vectors n± point respectively from the charge's 
advanced/retarded position to the position Rn on the 
sphere 523' . Trajectories are assumed to be bounded such 
that |xi:(tfc)| << i?, so that for each charge we have 
n+ ~ n_ = n. 

The retarded far-electric and far-magnetic fields of a 
charge qk at the space-time point (t, Rn) are piecewise- 
defined by the Licnard-Wiechert formulas |23j 



Tpret 



(t,n) 



gfe n X [(n - v^) x at] 



R 



{l-n-wkf 



(15) 



and 



Br*(i,n) = nxEr*(t,n). 



(16) 



In Eq. ([T5|) we have used the far-field-limit in which the 
light cone distance 



(17) 



rk{tk) = \yik{tk) - Rn\ 



is equal to R since |xfc(ifc)| « R. In Eq. pSj) 
Vfe = dxfc/dtfclf^ and a^. ee d^y^k/dtWi^ are respectively 
the charge's velocity and charge's acceleration at the re- 
tarded time ifc defined implicitly and piecewise by the 
retardation condition 



(18) 



tk^t- |xfe(ife) -i?n|. 



where | • | denotes Cartesian distance. Equation (|25|) de- 
fines tk as an implicit function of time t with a piecewise 
defined derivative 



dtk 
dt 



1 



(1 - n • Vk] 



(19) 



Using Eq. (jTS]) to evaluate the far-magnetic field ([T 
yields 



Br*(t,n) = -^x 



afc 



(n • afc)vfe 
^(l-n.vfc)2 + (l-n.vfc): 



(20) 



The trajectory Xfe(ife) is a function of t from Eq. (|18p so 
using the chain rule and Eq. (J19p twice we can re-write 
the far-magnetic-field ([20]) as 



(21) 



Br(i,n) = -^x4K(i,)] 



i? 



dt2' 



The far-electric field is a linear function of the far- 
magnetic field obtained using Eq. (|16|) and the transver- 
sality property n • E)J'^*(t, n) = of the far-electric field 
(HH), i.e., 



Er*(i,n) = -nxBr*(i,n). 



(22) 



In view of Eq. ((22)) . it suffices to study the vanishing 
of the retarded-far-magnetic-fields. We further assume a 
symmetry that the time-reversed orbit yields the same 
orbit rotated about an axis. For these reverse-rotate- 
symmetric orbits the vanishing of the retarded far-fields 
implies the vanishing of the advanced far-fields. 

From now on charge 1 is taken to be positive and equal 
to q while charge 2 is negative and equal to —q. The GAH 
along a bounded piecewise C^ orbit is then expressed 
almost everywhere by 



B 



Bl" 



B-* 



on d , , . 
x-(x,(.0 



R 



X2(t2)) = 0. 

(23) 

In the family of orbits with discontinuous velocities one 
can readily construct bounded orbits with vanishing far- 
fields; e.g. piecewise-constant-vclocity orbits with trajec- 
tories consisting of polygonal lines. These are bounded 
orbits with each acceleration vanishing piecewise, so that 
the radiation fields vanish. The question that needs to be 
answered is "Do we need these velocity discontinuities?" . 
Equation (j23p is a (linear) neutral differential delay 
equation with piecewise-linear continuous solutions de- 
fined on the intervals t G (tjr_j^,t^), with it € Z by 

Xi(ii)-X2(i2) =D,(n)+n/,(<,n) + (i-t-)-K.(n), (24) 

where the D^ (n) and Va (n) are arbitrary bounded func- 
tions and the /o-(i, n) are bounded and piecewise C^. It 
is possible to choose n • 0^(11)= and adjust D(j(n) in 
each interval to make the left-hand-side of Eq. ([M]) con- 
tinuous. 

Along a spatially bounded orbit, Eq. ([T51) is approxi- 
mated for large values of R by 



tk^t - R + n ■ Xk(tk 



(25) 



Notice that Eqs. (j25l) yield an implicit relation between 
ti and ^2, 



ii - <2 = n • (xi(ii) - X2(t2))- 



(26) 



Given the trajectories Xi(ti) and X2(i2), Eq. ([26l) and 
the implicit function theorem yield ii as a function of ^2 
and n. Define the influence interval of point (^2^x2(^2)) 



by the interval containing ii when n varies arbitrarily in 
Equation (|26)) . i.e., 



i2-|xi(tl)-X2(t2)| <ti <t2+|xi(ii)-X2(i2)|. (27) 

The time span (l?f)) is from the retarded lightcone- 
time of (^27X2(^2)) to the advanced lightcone-time of 
(t2,X2(t2)), as along the sewing chain illustrated in Fig. 
1. Notice that the future lightcone appeared naturally 
in the two-particle problem, even though we were deal- 
ing only with the retardation conditions (|25p. It follows 
from Eqs. ^ and ^ that 

/<,(i, n) = (ti - i2) -{t- t-)n ■ V„{n). (28) 

and we can therefore re-write Eq. (IMl) as 

xi(ii) - X2(i2) = D<,(n) + (ii - i2)n -{t- t-)n x L^(n). 

(29) 
where Lc,(n) = nxl/o-(n). Since linear growth in a con- 
stant direction is unbounded, the only globally C^ orbit 
must have £(j(n) = Vcr, and it follows from Eq. ([29]) 
with £cr(n) — that 



xi(ii) - X2(i2) = D„(n) + (ii - t2)n. 
The derivative of Eq. (l30l) respect to time yields 



vi 



V2 



(1-n-vi) (l-n-V2) 



Ki2n 



where 



K,o = 



1 



1 



(1-n-vi) (l-n-V2)' 



(30) 



(31) 



(32) 



Equation (j26l) allow us to move n in a cone with axis 
along xi(ti) — X2(i2) 7^ in a way that fixes ii and ^2 
while changing t with Eqs. ([25)1 . On the other hand, for 
fixed ii and ^2 the left-hand-side of Eq. ([3T|) spans a 
plane of the fixed vectors vi(ti) and V2(t2), so that Eq. 
([31]) can hold only if K12 = 0, which combined with Eqs. 
([3T1) and dSS]) yields 



Vl(tl) = V2(t2)- 



(33) 



Equation (|33|) defines globally-constant velocities along a 
fixed direction, which in turn implies unbounded motion 
unless vi = V2 = 0, as discussed in Ref. [12]. This 
impossibility follows if velocities are to be continuous. 

Nontrivial alternatives to this unsatisfactory conclu- 
sion necessitate the introduction of discontinuities by 
varying the direction of the piecewise- velocity-like term 
io-(n) ^ of Eq. (|29l) in each interval. The piecewise 
derivative of Eq. ([^ respect to time yields 



vi 



V2 



(1-n-vi) (l-n-V2) 



i4:i2n-nxL^(n). (34) 



Notice that K12 is still given by Eq. ([5^ and with 
nonzero L^ija) the right-hand-side of Eq. (IMl) forms a 



complete 3-dimensional basis to express any vector (in- 
side or outside the plane of Vi(ii) and V2(i2) )■ Equation 
(PS)) still allows one to move n in a cone with axis along 
xi(^i) — X2(i2) ^ in a way that fixes ti and ^2 while t 
changes with Eqs. (P5|) . By choosing K12 and a nonzero 
Lcr{n) for each t G {t'^_^,t'^) we can describe any vec- 
tor on the left-hand-side of Eq. (l34l) , so that there is no 
inconsistency. 

As an example, time-reversible orbits satisfying Eq. 
([34|) are piecewise-constant-vclocity orbits generated by 
having one velocity jump at a given time while the 
other velocity jumps either in the backward or forward 
lightcone-times symmetrically, as well as at every time 
in the forward and backward light-cones of a discontinu- 
ity time (the sewing chain illustrated in Fig. 1). These 
piecewise-linear polygonal orbits can be shown to satisfy 
Eq. (1^ by direct substitution and use of Eq. ([TO]). 
In the following we show that Eq. ((29|) and the implicit 
function theorem yield a consistent piecewise-defined tra- 
jectory Xi(ii) from a given piecewise-defined trajectory 

X2(t2). 

Notice that for given continuous and piecewise C^ 
X2(i2), Do-(n) and Ldn), in general Eq. ([^ deter- 
mines only a function Xi (ti , n) of the two variables (ii , n) 
through 

Xi(ti,n)=X2(t2)+D^(n) + (ti-t2)n-(f-t-)nxL^(n). 

(35) 
The implicit function theorem further determines ^2 and 
t as functions of ii and n from Eqs. ((25|) . (|26p and ([35]). 
For the implicit function theorem to yield a consistent 
trajectory, we must satisfy the consistency requirement 
that xi(ti, n) determined by Eq. ((35|) is a function of ii 
only, i.e.. 



9xi(^i,n) 
dn 



0. 



(36) 



Condition ([M|) applied to the right-hand-side of Eq. 
(|35p is the extra condition determining a consistent tra- 
jectory. Since condition (pG]) must hold for all values of 
ii in each piecewise interval of the orbit, we must also 
have inside each piecewise interval that 



9^X1 (ti,n) _ d 9xi(ti,n) 



dtidn dn dti 

which can be expressed as 



-)=0, (37) 



d ,, 9<2(ii,n) at(ti,n) , ^ m n 

^[(V2 - n)^^— + n - ^^„ X LM)\ = 0- 

The general solution to Eq. psp involves an arbitrary 
piecewise-defined function ^0.(^1)1 i-e., 

dt2{ti,n) dt{ti,n) 
nH ^r ^(v2 - n) '— — -n x i^(n) = A„{ti). 



A symmetric condition follows by exchanging indices 1 
and 2 in Ea. (l38l) . introducing an arbitrary i?o-(i2) and 
changing the sign of L^ (n) , yielding 

dti(t2,n) , . dt(t2,n) ^ , ^ „ , n 

n+ 'L ^ (vi - n) + ^'' ' n x i,(n) = ^,(^2). 



dh 



dto 



(40) 

The partial derivatives in Eqs. ([39l) and (|40l) can be 
evaluated using Eqs. (|25|l . yielding 



= Xi2n-nxL,(n)(41) 
= Xi2n-nxL,(n)(42) 



(1-n-vi) (l-n-V2) 
vi(ii) B.(t2) 



ail 



ati 



(39) 



(l-n-vi) (l-n-V2) 

where again K12 is defined by Eq. ([32|) . 

Notice that Eqs. (gl]) and (j42l) with ^^.(ii) = vi(ii) 
and Bc{t2) = V2(t2) bring back Eq. ([M)) as the single 
necessary condition to construct a consistent piecewise 
C^ continuous trajectory xi(ii) from a given piecewise 
C^ continuous trajectory X2(t2) by the implicit function 
theorem. We stress that velocity discontinuities are ab- 
solutely necessary. The discontinuities introduced by the 
nonzero L„ (n) in each piecewise interval play an essential 
role in the solution as discussed below Eq. (p4)) . Trying 
to solve either Eq. (jH]) or Eq. (|42|) with i<j(n) = 
stumbles with the former obstruction that rotation of n 
with fixed ti and ^2 places n either outside the plane 
of Acr{ti) and V2(i2) or outside the plane of Ba{t2) and 
vi(ii). The nonzero L„{r\) provides the third linearly in- 
dependent direction forming a complete basis to express 
nin Eqs. (gl]) and ([42]) . 

It can be seen that piecewise-constant-velocity polyg- 
onal orbits can be constructed using the suitable Lo.(n) 
defined by Eq. (j34p . after which the implicit function 
theorem constructs consistent continuous piecewise C^ 
trajectories. It would be desirable to find bounded min- 
imizer orbits of the variational method [9] satisfying the 
vanishing far-field conditions ([1T|) and (g^. as first con- 
jectured for the orbits studied in Ref. [18]. The justi- 
fication to generalize to trajectories with discontinuous 
derivatives is to include short-range bounded GAH or- 
bits in the family of physically possible orbits. This was 
used in Ref. [18| to predict spectroscopic lines of hydro- 
gen within a few percent agreement with the predictions 
of quantum mechanics. 



DISCUSSION AND CONCLUSION 

The fact that accelerations are discontinuous is ex- 
pected because the Wheeler-Feynman equations of mo- 
tion are explicitly neutral for the accelerations. Con- 
sequently, it could seem that a theory of piecewise- 
restricted Wheeler-Feynman equations of motion should 
have only acceleration discontinuities, a fact that already 



8 



introduces discontinuous fields and demands a general- 
ization of electrodynamics. We have seen that general- 
izing to trajectories with discontinuous accelerations is 
not sufficient to include bounded two-body orbits with 
vanishing far-fields. Our analysis starting from the vari- 
ational method as the fundamental principle has shown 
that, in general, the velocities are also expected to be 
discontinuous at the same "generalized breaking points" 
along the minimizer orbits. Our analysis, using the vari- 
ational method as a boundary- value problem, shows that 
the most general solution of the Wheeler-Feynman neu- 
tral differential delay equations has discontinuous accel- 
erations and velocities. 

The form of Eq. (|12p is reminiscent of the energy oper- 
ator used formally in quantum mechanics. Since the evo- 
lution parameter is arbitrary, the same parameter can be 
used in a Lorentz-transformed frame, such that the mo- 
mentum currents transform by like a four-vector. The 
existence of four components that must be continuous 
and given by the partial derivatives of a scalar invari- 
ant is again analogous to the quantum Dirac equation 
and suggests a property analogous to spin for the point 
charges. It is remarkable that Wheeler-Feynman elec- 
trodynamics completed with a finite action endows the 
point charges with a spin-like property. The existence 
of a spinorial four-component momentum current (Eqs. 
PTjl and (|12l) ) is due to the parametrization-invariance 
symmetry of the electromagnetic variational method [^ . 
Otherwise a generic action with delayed interaction has 
only three momentum currents continuous at breaking 
points of minimizer orbits. 

In Ref. [29[ only globally C^ solutions were sought 
for the seemingly non-neutral one dimensional motion, 
so that piecewise-defined solutions with discontinuous 
velocities awaited study. In considering this, it is im- 
portant that the electromagnetic-action-functional ([T]) 
of Ref. [9| yields a neutral-boundary-value-variational- 
method, as opposed to the non-neutral variational meth- 
ods of classical mechanics. 

The variational principles of classical mechanics are 
two-point boundary value problems that are equivalent 
to an initial-value problem with initial velocity chosen to 
hit the final trajectory point. Moreover, in the classi- 
cal problem there is no issue of velocity continuity, be- 
cause the "history" for a finite-dimensional ODE is a 
point. On the contrary, for the electromagnetic varia- 
tional method, choosing the "initial velocity" to shoot 
the final point either requires a velocity discontinuous 
with the past boundary history or the trajectory arrives 
at the final point with the wrong velocity. The gener- 
alization to discontinuous velocities extends the solvabil- 
ity of the electromagnetic-boundary- value-problem to a 
larger class of boundary value data, which is the second 
advantage of our extension of Wheeler-Feynman electro- 
dynamics. 

Electromagnetism was originally formulated with the 



integral laws of Ampere, Gauss and Faraday, and only 
much later differential equations holding everywhere were 
introduced by Maxwell. The requirement of a second 
derivative existing everywhere is actually not needed for 
particle dynamics, where one is concerned only with the 
integral of the force along the trajectory. The variational 
method is a step back from Maxwell's equations in the 
sense of weak solutions. Replacing Maxwell's equations 
by the vanishing of the Frechet derivative (jG]) along a 
continuous trajectory with boundaries in future and past 
yields solutions defined only on bounded time-intervals. 
From these segments of orbits one can construct the fields 
as derived quantities [8|, l22l | . Since fields constructed in 
this fashion involve a retarded and a advanced position, 
before defining fields everywhere in space we need to ex- 
tend to a global trajectory. Extension is possible using 
conditions (fTTj) and (fT2)) and in general involves velocity 
and acceleration discontinuities. It is then tempting to 
translate our generalized electrodynamic quantities into 
the concepts of Maxwell's electrodynamics, but it must 
be done carefully; The derivation of the differential form 
of Maxwell's equations given in [8| holds only in regions 
where the extended fields are C^ . Consequently, Poynt- 
ing's theorem in integral form holds in regions where the 
fields are C"^ . Nevertheless, quantities like energy fiux 
and field energy in a volume have a meaning even for 
discontinuous fields, so that a statistical interpretation 
could be sought in the case of discontinuous fields; For 
example, the generalized flux of the Foynting vector is 
an integral that can be evaluated even when fields are 
undefined on sets of measure zero. The Foynting vector 
P = ExB evaluated with Eqs. ^ and ^ at (t, Rn) 
becomes 



1 



{IE- 



adv\2 



IE' 



ret\2 



}n 



(43) 



where single bars denote the Euclidean modulus. 

We stress that a condition of non-radiation is weaker 
than our GAH, as follows; The GAH implies the vanish- 
ing of the flux (021) because \W*\'^ = jE'^'^^'p = almost 
everywhere, while the converse is not true, i.e., the van- 
ishing of the flux integral alone does not iniply the GAH. 
For example, the circular orbits of Refs. [13|, |lj| do not 
satisfy the GAH, and even though these orbits do not 
radiate on average, circular orbits have non- vanishing 
far-fields to disturb a "third charge of the variational 
method" (i.e., are not short-range). In Ref. IJ] a model 
for the neutron was attempted, and even if it had not 
failed for other reasons, it would yield a neutron with 
far-fields. As regards non-vanishing far-fields, a first at- 
tempt to overcame the GAH-deficiency of circular orbits 
14 1 was the perturbation theory of Ref. [18| that added 
high-frequency modes of the tangent dynamics to enforce 
Eq. ([23l) at the frequency of the circular orbit. 

We have shown that our variational method yields a 
dynamical system even for limiting orbits with discontin- 



uous velocities. For example, along piecewise-constant- 
velocity polygonal orbits the variational equations of mo- 
tion would be applied as follows; On discontinuity corner 
points with no acceleration defined, one enforces conti- 
nuity of momentum only, Eq. (|lip. At other points, 
wherever accelerations are defined, one uses the usual 
Wheeler-Feynman equations of motion. (Notice that 
piecewise-constant-velocity polygonal orbits have vanish- 
ing far-fields but obviously do not satisfy the equations 
of motion, unless charges are far apart). 

We have demonstrated five different and important 
reasons to study orbits with discontinuous derivatives: 
(i) inclusion of bounded GAH orbits as short-range or- 
bits; (ii) compatibility with the conservation of Noether 's 
momentum; (iii) compatibility with the neutrality of the 
equations of motion of the Wheeler-Feynman electrody- 
namics; (iv) the fact that the variational method is natu- 
ral in a space completed to contain orbits with discontin- 
uous velocities; and (v) inclusion of limits of sequences 
of non-radiating orbits. 

The physical need for trajectories with discontinuous 
velocities is justified as limiting orbits defined by Cauchy 
sequences of bounded orbits which must develop kinks 
in the short-range limit (i.e., the GAH). In Ref. 18 1 
the GAH-deficiency of circular orbits was removed with 
a perturbative Fourier series that solved Eq. (p3| at the 
first harmonic frequency only. As we have shown here, 
the perturbative series of Ref. 18| should converge to 
an orbit with discontinuous velocities. The short-range 
condition of Ref. |18| predicted orbits and spectral lines 
in the atomic magnitude with a surprising precision, so 
that we can claim agreement with experiment and quan- 
tum mechanics. Piecewise-defined minimizers have also 
been used successfully to explain double-slit diffraction 
in Ref. [I^. 

From our generalized electrodynamics with discontin- 
uous derivatives it should be possible to derive a general- 
ized electrodynamics with delayed-only interactions and 
self-interaction, using the generalized absorber hypothe- 
sis (GAH) in close analogy with the derivation of Wheeler 
and Feynman |8i . Since we expect solutions with velocity 
discontinuities, Taylor expansions of deviating arguments 
should be avoided or piecewise-restricted. It is known 
that many-component delay differential equations behave 
like neutral differential delay equations when some solu- 
tion components are discontinuous at breaking points, in 
the sense that the discontinuous derivatives never smooth 
out [3l[. Therefore, the third derivative should be gen- 
eralized by restricting it to a left-derivative and a right- 
derivative at breaking points. Also, the generalized ab- 
sorber hypothesis with discontinuous fields no longer im- 
plies the vanishing of the difference of retarded and ad- 
vanced universal fields everywhere, as used by Wheeler 
and Feynman [8|, |32|. Our Eq. (|34|) and its advanced 
version give the corresponding weaker generalization to 
this former stronger condition of a vanishing difference 



of retarded and advanced fields everywhere. We spec- 
ulate that Eq. (|34p should be the starting point for a 
generalized theory of self-interaction free of the perva- 
sive runaways of the two-body problem with the usual 
self-interaction. 

Last, we speculate that unbounded scattering orbits are 
different from the bounded orbits studied here. Along un- 
bounded orbits Eq. ([M)) contains the extra secular term 
with a constant V^(n) ^ 0. The dependence on bound- 
ary segments and time separation must be investigated 
for scattering trajectories with discontinuous velocities 
and accelerations at the boundaries; For example, even 
if the history segment (0~, O"*") of Figure 1 is assumed 
C°°, the forward sewing chain of 0+ places a breaking 
point /a in the history segment (£^, L'^) of particle 1. 
Unless histories are very special so that derivatives are 
continuous at O"*", in general the history (L~, L+) should 
involve a discontinuous derivative at point /s. Scattering 
trajectories are likely to have future continuations involv- 
ing stiffer jumps at later times, so that particles collide 
with laboratory boundaries, which can be regarded as a 
generalized type of radiative loss. 



Figure Captions 



Figure 1: Illustrated in red is the initial point Oa of 
trajectory 1 plus the segment of trajectory 2 inside the 
lightcone of Oa, i-e., from point 0~ to point O"*" and the 
endpoint Lb of trajectory 2 plus the segment of trajec- 
tory of particle 1 inside the lightcone of Lb, i.e., from 
point L~ to point L^ . The trajectory of particle 1 of the 
variational method goes from Oa to L~ (blue line) while 
the trajectory of particle 2 goes from 0+ to Lb (green 
line). The first breaking point is point O"*" which gener- 
ates a forward sewing chain of breaking points /i, /2,/3 
while endpoint L~ is a breaking point generating a back- 
ward sewing chain of breaking points 6i, 62, ^3- 
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